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SUMMARY 


Owing to the quantity of material to be presented, this 
report is divided into several volumes. Volume 1 introduced 
the series, described the analytical procedures in detail, and 
presented and discussed the results for Shot 10. Volume 2 
presented and discussed the results for Shot 9. This volume 
presents and discusses the results for Shot 8. A subsequent 
volume will present the results for Shot 11, and compare the 
results of the four experiments. The method of analysis is 
common to all four experiments and is described in detail in 
Volume 1 only. 

So that the results from the four experiments may be 
easily compared, they have been scaled to remove the effects 
of varying atmospheric conditions. (Results are scaled to a 
1 kg charge weight and a standard atmosphere of dry air at 
15°C at sea level.) For the most part, only scaled results 
are presented. Exceptions include some derived pressure-time 
histories, which are compared to actual gauge measurements 
made in the experiment. 

Results are presented in SI units, even though the experi- 
ments were originally laid out in British units. Only distance 
and time measurements are affected, however, as velocity, 
density, and pressure results are presented as dimensionless 
ratios. A distance units conversion scale is included 
to convert between SI units (meters scaled to a 1 kg charge) 
and British units (feet scaled to a 1 lb charge), plus a time 


scale factor and scale factors to convert pressure ratios to 
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both British and SI pressure units. Scale factors which may 
be used to compute the distance and time values actually 
observed under the ambient conditions of each shot are also 
provided. Dimensional pressure units are used for the results 


presented at gauge locations. 
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CHAPTER 1, SHOT 8 ANALYSIS 


1.1 Introduction 


This is the third volume in a series which presents the 
particle trajectory analysis results from four experiments 
(Dipole West Shots 8, 9, 10 and 11) carried out to obtain 
information on the interaction of spherical blast waves with 
real and ideal reflecting surfaces. A general description of 
the project can be found in Volume 1. The results presented 

in this volume are for Shot 8. 

In each experiment, photogrammetrical studies were made 
of the shock fronts (refractive image analysis, RIA), and of 
the motions of smoke puff particle tracers (particle trajectory 
analysis, PTA). The refractive image analysis results were 
reported by Dewey et al. (1975) and results of the particle 
trajectory analysis are presented in this report. The method 
of particle trajectory analysis, common to all four shots is 


described in detail in Volume 1 only. 
1.2 Description of Shot 8 


Dipole West Shot 8 was fired on September 17th, 1973 by 
the Ballistics Research Laboratories at the Defence Research 
Establishment Suffield, in Alberta, Canada. Two 1080 lb (490 kg) 
spheres of Pentolite were detonated simultaneously, to within 


5 microseconds, at nominal charge heights of 25 and 75ft 


(7.6 and 22.9m) over a relatively smooth ground surface. 


Particle trajectory data were gathered by photographing 
the movement of smoke puffs formed in a vertical plane running 
out from ground zero at 6.7° south of west. A WF5 camera 
operating at about 3500 frames per second was positioned 
30ft (9.2m) above ground level at a position 600ft (183m) 
due south of ground zero (GZ), the point on the ground 
vertically beneath the charges. 

Table 1 gives the field survey data for the event, and 
Figure 1 shows a plan view of the layout. The dashed line 
represents the approximate line of sight of the WF5 camera. 
Figure 2 shows the field of view of this camera. 

The smoke puff grid was made up of 20 columns of 12 puffs 
each, hung vertically on strings. The vertical spacing of 
puffs was 5ft, beginning 3ft above ground level and ending 
at a height of 58ft. The horizontal spacing of the columns 
of puffs was 10, 7 or 5ft, depending on the distance from 
ground zero, beginning at about 25ft and ending at about 140ft 
from GZ. Of the possible 240 smoke puffs, 237 detonated 
successfully. A good film record was obtained. 

This report describes the analysis of the smoke puff 
data collected for Shot 8, and presents and discusses some 


of the results of that analysis. 


1.3 Camera calibration and data reduction 


The calculated camera position coordinates and orien- 


tation angles for Shot 8 are presented in Table 2, together 


with the positions of photomarkers transformed from one frame 
of the film just before detonation to an object plane defined 
as passing thiough ground zero and being normal to the camera 
orientation axis. The differences ("shifts") between the 
object plane positions of the transformed calibration points 
and their positions computed from the field survey data are 
given in Table 2. The object plane positions of the calibra- 
tion points computed in these two ways are also shown in 
Figure 3. 

The camera calibration procedure, described in detail 
in Volume 1, ensured that selected photomarker images (Pl to 
P5) transformed to the object plane in a way which matched 
them exactly to the positions computed using the survey data. 


These reference photomarkers for Shot 8 are indicated in Figure 3 


using large circles: namely, Pl = Wl, P2 = W3, and P3 300W1. 


The separation distance between P4 = P3 300W1l and P5 300wW2 
was also used as a calibration parameter. Photomarker VPI1A 
was missing on Shot 8. 

The image positions of two reference photomarkers 
(VP3B and 300W2) and all smoke puffs were measured frame~-by- 
frame over a time interval corresponding to the approximate 
duration of the positive phase of the blast waves (film frames 
9 to 375), and were transformed to distances in the object 
plane by matching the reference marker positions to their 


positions transformed from the calibration frame. These data 


were again transformed from the object plane to the smoke puff 


det 


a 


plane which was assumed to pass through "corrected" ground zero; 
to be vertical, and to run 6.7° south of west from GZ. 

The x-y coordinate system in the smoke puff plane was 
the same for Shot 8 as for Shots 9 and 10 except that the 
corrected value for ground zero was displaced 0.8ft from the 
surveyed ground zero, in a direction approximately 59° north 
of east. The corrected ground zero was defined to have the 
same elevation as the surveyed ground zero, but was located 
directly under the midway point between the two charge centers. 
As for Shot 9 and 10, all data in the output plane are plotted 
with the x coordinate reflected, i.e. with positive values of x 
to the right hand side, as if the smoke grid had run to the 
right of the charges rather than to the left as seen in the 
film images. 

A time was assigned to each film frame using the 1 ms 
timing marks placed on the film during its exposure. The 
film timing method was described in Volume 1, and the complete 
set of film timing data used for Shot 8 is provided in Table 3. 

Figure 4 shows the positions of the detonated smoke puffs 
at a time prior to the detonation of the two charges. These 
positions are in the plane of the charges and the smoke puff 
grid, as described above. The smoke puff plane was not exactly 
parallel to the camera image and object planes (Figures 2 and 3), 
and various geometrical corrections were applied to make the 
transformation between them. The puffs enclosed in parenthesis 
were not visible in the earlier film frames because they were 


concealed by photomarkers, but were seen later. Charge 
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positions in the figures are plotted as if they were 
positioned exactly above the corrected ground zero origin. 


The data shown in Figure 4 have not been scaled. 


1.4 Data scaling and trajectory fitting 


The position-time histories of individual smoke puffs were 
extracted from the frame-by-frame positions of the smoke puff 
grid, and then scaled to standard atmospheric conditions and 
charge weight. A change to SI units was made at this point in 

the analysis. The resulting trajectories were edited, and 
then smoothed by fitting polynomial functions. 

Particle trajectory data were scaled by dividing all 
distances by Sachs scaling factor S = V(WP OY 7 (WoP) and 
multiplying all times by the factor C/(C5S), where C is the 
ambient sound speed computed for Shot 8. Data used to compute 
C and S, and define the scaled event, are listed below with the 


computed values of C and S. 


Ambient temperature, Pa= WOT 72" oC (67.5 °F) 
Ambient pressure, P = 93.22 kPa (13.521 PSL) 
Relative humidity, RH = 31.0% 

Computed vapour pressure, VP = 0.71 kPa (5.3 mm Hg) 
Ambient sound speed, C = 343.635m/s (dilz7 2tys) 


Charge weight, W = 489.9 kg (1080 lbs) 


Sachs scaling factor, S = 8.1051 

Standard charge weight, Wo = 1.0 kg {2.2 lbs) 

Standard pressure, Py = 101.325 kPa (14.7 PST) 

Standard temperature, a 5S le (59 °F) 

Standard sound speed, C_ = 340.292 m/s (1116 ft/s) 
(dry air) 2 


The results presented in this report, therefore apply to 
a scaled event which is the detonation of two 1 kg charges in 
a standard atmosphere. The scaled heights of burst for Shot 8 
were 0.919 m and 2.795 m, and the scaled charge separation 
divided by two, was 0.938 m. 

Figure 5 shows the scaled particle trajectory data for 
Shot 8 in the smoke puff plane with positions measured horizontally 
and vertically from corrected ground zero. Approximately 
26509 puff positions are represented. As represented, the raw 
trajectory data have not been smoothed. 

The raw particle trajectory data were edited to remove 
obvious data processing errors, such as a Single point widely 
displaced from its trajectory for one or two frames. The 
trajectory of each puff in turn was then smoothed by least 
squares fitting simple polynomial expressions separately to 
both the x and y coordinate data, these being discrete functions 
of frame time. The adequacy of each fit was determined by 
examining on the same graphical output, plots of both the raw 
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trajectory data and the fitted curve. For Shot 8 this meant 
examining and adjusting 474 such plots, at least two or three 
times each. 


For a given puff, 


the first step in fitting the raw 


trajectory data was to set the time of arrival of the shock 


front first hitting the puff. The data at subsequent times were 


fitted with polynomial functions, as described in Volume l, 


paragraph 2.5. The first derivatives of the fitted functions 
were also calculated at a series of times for use in later 


calculations of particle velocity. 


1.5 Regionalization and shock strength calculations 


Five regions were defined in the smoke puff plane on 


the basis of the shock front which first struck the puffs in 


a particular region. 


These are shown in Figure 6. The regions 


were bounded by the triple point trajectories measured using 


refractive image analysis (Dewey et al., 1975). Regions 
1 and 2 are those in which the smoke puffs were first hit by 
a spherical primary shock front, and regions 3, 4, and 5 are 


those in which the puffs were first hit by a Mach stem. 


In each of the five regions, the shock trajectory data 


obtained from the first movement of the smoke puffs were fitted 


to a function of the form 


: r(t) = A + Bt +C log (1 + t) + D/ log (1 + t), 


where r is the shock radius, t is the time after detonation, 


and A, B,C andD are the fitted coefficients. 


The shock 
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velocities were calculated by differentiating this function. 
The peak particle velocity, ae , peak density, Der and peak 
hydrostatic overpressure, Por as functions of shock radius in 
each of the five regions, were calculated from the shock 

q velocity using extensions of the Rankine-Hugoniot equation. 


Details of the shock radius calculations etc. are described 


in Volume 1, paragraph 2.6. 
1.6 Particle velocity calculations 


Particle velocities were computed using the methods 


described in Volume 1, paragraph 2.7. 
1.7 Density and hydrostatic overpressure calculations 


Densities and hydrostatic overpressures in the smoke puff 
plane were calculated by the method described in Volume l, 
| paragraph 2.8. Results in both cases represent average 


values over cells defined by four adjacent smoke puffs. 
1.8 Surface representation 


Surfaces were fitted to the times of shock front arrival 
and to the fields of particle velocity, density and hydrostatic 
overpressure at a sequence of times. All data were interpolated 
onto a common regular Euleurian grid. Fields of dynamic pressure 
were computed from surface-interpolated particle velocity and 
density results. Contour plots were generated for all 


surfaces at selected times, and time histories computed at 
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several fixed locations. The methods used were identical to 


those described for Shot 10 in Volume 1, Chapter 3. 
1.9 Pressure and total pressure time-histories 


To permit a direct comparison between results obtained 
from the particle trajectory analysis and measurements made 
using side-on and face-on pressure transducers, the hydro- 
static and total overpressure time-histories were calculated 
at those locations coincident with guage positions within the 
smoke puff grid. Dynamic pressures and hydrostatic over- 
pressures obtained from the particle trajectory analysis were 
used to compute the total pressures after application of a 
compressibility correction. This correction is a function of 
the local Mach number and its form depends on whether the Mach 
number was greater or less than unity. The time varying hydro- 
static and total overpressure impulses, determined by integrating 
the pressure time histories, were also calculated and compared 
with similar integrations of the electronic transducer data. 

The methods used to calculate the total pressures and the 
impulses are described in detail in the addendum to Volumes 
1 and 2, which is incorporated in this volume. In cases where 
the leading edge of a time-history curve was rounded, integra- 
tion of impulse was done using data interpolated linearly 
between the peak parameter value determined at the time of 


arrival, and a point on the time-history curve subsequent to 


the time of arrival. The second of these two points was chosen 
17 


in a manner which ensured a minimum difference in slope between 


the interpolated and computed sections of the time-history data. 
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CHAPTER 2. SHOT 8 RESULTS 
2.1 Times of shock front arrival 


The measured initial puff positions, the times of first 
shock arrival, and the peak particle velocities obtained by 
differentiating the functions fitted to the particle 
trajectories are presented in Table 4. Puff position is 
given relative to corrected ground zero as origin, with hori- 
zontal and vertical axes. Puff position and the time of 
arrival of the first shock are given both as observed and 
scaled. Particle velocities listed are derivatives of the 
fitted puff trajectories at the times of shock arrival, and 
are expressed in Mach units. Expressed this way, the particle 

velocities are the same scaled as unscaled. Also listed are 
the initial radial puff positions (scaled) and region codes. 
Shock front data determined from the first movement of 
the smoke puffs, i.e. calculated from the tir --of-arrival 
data in Table 4, are listed in Tables 5.1 - 5.5. Each table 
corresponds to one of the 5 regions used. Listed are the 
observed and fitted unscaled shock trajectory data, the 
scaled fitted shock trajectory data, and the computed shock 
velocities and peak parameters associated with shock strength: 


peak hydrostatic overpressure in atmosphere and kilopascals, peak 
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particle velocities in Mach units, and peak density ratios. 
Given as ratios, these peak parameters are the same scaled 
as unscaled. Pressure given in kilopascais in the tables 
refers to the unscaled (observed) case only. 

The shock front radius versus time data derived using 
particle trajectory analysis (PTA) are also shown in Figures 
7.1 - 7.3 for the two primary fronts, the two Mach stems at 
the interaction plane, and the ground Mach stem, respectively. 
They are compared to corresponding data derived from refractive 
image analysis (RIA) reported by Dewey et al. (1975). The 
refractive image analysis results were obtained using photo- 
grammetry against a striped canvas backdrop and they describe 
the shock as it travelled in a direction almost diametrically 


opposite to the direction of the smoke puff grid. 


2.2 Shock strengths 


Peak particle velocities calculated from shock front 
velocities are shown inFigures 8.1 - 8.3 for the primary 
fronts, interaction Mach stems, and the ground Mach stem. 

This method of determining peak particle velocities has been 
labelled method 1, and the data plotted correspond to those 
listed in Tables 5.1 - 5.5. The results in the figures are 
compared with those previously obtained using refractive 

image analysis (RIA). In the case of the primary shock fronts, 


results are also compared to those of Brode (1957) for TNT. 
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In Volume 1 other methods of determining shock strengths 
in the various regions were described. It was demonstrated 
that method 1 was clearly the most accurate, and in the present 
volume shock strengths calculated using methods 2 and 3 are 
not reported. For this reason Figures 9, 10 and 11 and 


Table 6 do not appear in this volume. 
2.3 Particle velocity fields 


The calculated particle velocities in the plane of the 
smoke puffs are shown as vectors in Figures 12.1 through 12.9, 
for various times after the detonation. All times and 
positions are scaled to a l kg charge in a standard atmosphere. 
The particle velocity vectors represent the derivatives of 
the smoothed particle trajectories, and their magnitudes may 
be judged using the standard vector shown on each figure. All 
velocities are measured in Mach units, relative to the standard 
sound speed. Puffs not yet struck by a shock wave are 
represented by small circles (zero velocity). 

Numerical data corresponding to Figures 12.1 - 12.9 are 
listed in Tables 7.1 through 7.12, along with scaled radial 
positions of the puffs, and region codes as defined in Figure 6. 
Conversion factors are given at the foot of each table, which 
may be used to convert the scaled data in the tables and figures 


back to their original unscaled values. 
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2.4 Density and hydrostatic overpressure fields 


Calculated average relative densities throughout the 
smoke puff plane are depicted graphically in Figures 13.1 - 13.4, 
for various times after the detonation. All time values are 
scaled. Cell positions are scaled and are given relative to 
the corrected ground zero as origin with horizontal and 
vertical axes. The calculated densities may be judged using 
the density shading scale shown on each figure. Density is 
given as a ratio, relative to ambient density. Cells not yet 
struck by a shock wave and cells in which the density has 
dropped to a value less than ambient density are shown blank. 

Corresponding numerical data are listed in Tables 8.1 - 
8.9 along with radial cell positions computed according to 
the regions defined previously. Numerical data describing the 
fields of hydrostatic overpressure are similarly listed in 
Tables 9.1 - 9.9. The pressure results for a given cell were 
obtained by multiplying the density results for that cell by 
a factor determined by the strength of the shock which first 
traversed the cell and which then remained constant, i.e. by 


assuming isentropic flow after the first shock. 
245 Times-of-arrival surface 


Figure 14 shows a perspective view of the surface fitted 
to the original smoke puff positions and the observed times of 


first shock front arrival, i.e., to the data listed in Table 4. 


ae 


The grid mesh size is 0.1 by 0.1 meters (scaled), about 


2.5 feet square (unscaled), or about % that of the original : 
smoke puff grid. The charge positions are indicated on the 
vertical distance axis. 


The times-of-arrival surface is smooth enough to permit 


contouring, the contours in this case (isochrones) representing 
shock front shapes at different times, as shown in Figure 15, 
but the surface is not smooth enough to permit the calculation 
of gradient vectors which could be used to compute shock 
velocity vectors and shock strengths over the new grid. 

Two attempts were made to obtain contours of shock 
strength. In the first, the times-of-arrival surface was 
smoothed by least-squares fitting low-order, one-dimensional 
polynomial functions to the time-of-arrival data along each 
grid row and column separately, and computing the derivatives 
of the fitted functions to obtain the associated components 
of the surface gradient vectors. Shock velocity vectors were 
obtained from the time-of-arrival gradients, and from these 
peak particle velocities were computed. The peak particle 
velocity (shock strength) surface is shown in Figure 16. 

The contours of this surface (not shown) did not exhibit any 
discontinuities across the boundaries of the shock front 
regions, as they would if surfaces had been fitted to the time 
of arrival in each region separately. 

The results of a second method used to compute shock 


strength contours are shown in Figure 17. These were obtained 
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by interpolating shock radius at each value of peak particle 
velocity shown, for each shock front region in turn, using 
the peak particle velocity versus radius curves shown in 
Figure 8. Arcs of circles with these radii, centered on the 
appropriate points along the vertical charge axis, were then 
drawn in the regions to represent shock strength contours. 
These peak value contours are discontinuous across triple 
point locii and other region boundaries. As a result, some 
horizontal lines are crossed twice by the same contour or, in 
other words, indentical shock strengths can be found at two 
locations the same vertical distance from a reflecting surface, 


but at different radial distances from the vertical charge axis. 


2.6 Field surface contours 


Contours of equal particle velocity, density, hydro- 
static overpressure, and dynamic pressure in the blast waves 
were determined for a series of times, using surfaces fitted 
to the various measured data fields at those times. Sample 
results are shown in Figures 18 through 21 at scaled times 
of 2.5, 4.0 and 9.0 ms. The shock fronts shown in these 
figures are obtained from the time-of-arrival surface (as were 
those in Figure 15). Field contours such as those shown can 
be drawn for any scaled time between 0.5 ms and 13.4 ms. 

It should be re-stated that all of these results were 


obtained from the photography of the smoke puffs only and do 


not rely on the results obtained using the refractive image 


analysis (Dewey et al., 1975). 
2.7 Time histories 


By mapping the physical properties of the blast waves 
at short time intervals its was possible to determine the time 
histories of these properites at any selected fixed position 
within the smoke puff grid. This was done at 15 fixed 
locations, three in the primary region of the lower charge and 
four in each of the three Mach stem regions, as shown in 
Figure 22. At each distance from the vertical axis through 
the charges in the Mach stem regions, each of the time 
history stations is approximately the same distance from either 
the interaction plane or the ground plane. (Particle velocity 
time histories could be interpolated closest to the grid edges 
because these were measured at puff locations, whereas the 
density and pressure data were measured at cell centers). 

Time histories of particle velocity, density, hydrostatic 
Overpressure and dynamic pressure at these locations are given 
in Figures 23 to 26. Time-histories of these physical properties 
of the blast wave can be provided at any other location within 
the smoke puff grid, on request. 

The vertical line which forms the leading edge of a time- 
history plot represents the interpolated time-of-arrival of 
the first shock at the given location, and the height of this line 
represents the peak parameter value derived from the shock velocity 


analysis. 
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The dynamic pressures plotted in Figure 26 are maximum 
values, computed using both the x and y component of 
particle velocity. Similar plots were made of the horizontal 
components of dynamic pressure, but the differences were not 
Significant since the y components of particle velocity at 
these locations were small. Other locations could have been 
chosen at which the y components would not have been insignifi- 
cant. 

Time histories for hydrostatic overpressure and total 
pressure are also plotted in Figures 27.1 to 27.4 for stations 
at the nominal positions of field-mounted pressure gauges on 
the "60 foot gun barrel". The gauges on this gun barrel were 
mounted at nominal elevations of 10, 20, 25, 40, 47, 50 and 
53 feet. The time histories at these locations are compared 
to the gauge measurements (Keefer and Reisler, 1975). The 
total pressures were calculated in the manner described in 
the addendum to Volumes 1 and 2 of this report. The variation 
with time of the integrated pressure (pressure impulse) is 
also shown in these figures, compared with similar integrals 


of the gauge data (Keefer and Reisler, 1975). 
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CHAPTER 3, DISCUSSION 


3.1 Particle trajectory analysis, Shot 8 


The methods used to analyze the smoke puff trajectories 
on Shot 8 were identical to those used for Shots 10 and 9 
and described in detail in Volumes 1 and 2 of this report. 
However, the results for Shot 10 clearly indicated the 
superiority of one of several methods of analyzing shock 
strength, and only the results of this method were reported 


for Shots 9 and 8. 


3.2 Primary shock strength of upper charge 


The refractive image analysis of the shock fronts 
described by Dewey et al., 1975 did not prove any infor- 
mation about the primary spherical shocks from the upper 
charges, and it was assumed that these charges had deton- 
ated satisfactorily. This assumptionwas validated for 
Shots 10 and 9 by the analysis of the particle trajectory 
time-of-arrival measurements. In Figure 7.1 the primary 
shock radii are plotted versus time for the upper and lower 
charges, for Shot 8, together with the results obtained for 
the lower charge by the refractive image analysis. All 


three curves appear to be identical. Unfortunately, because 


| 
| 
| 


of the limited range of the data obtained for the primary 
shock from the upper charge it was not possible to calculate 
accurately the variation of the shock strength with distance. 


However, these results for the lower charge, in Figure 8.1, 


show an identical shock strength variation with distance to that 
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obtained from the refractive image analysis and one which is 


very similar to Brode's (1957) calculations for TNT. 


3.3 Comparison of Mach shocks over different surfaces 


The refractive image analysis of Shot 8 (Dewey et al., 
1975) showed what appeared to be a small but significant 
difference between the strengths of the Mach shocks over the 
smooth ground and beneath the interaction plane between the 
two charges. The results of the particle trajectory analysis 
given in Figures 8.2 and 8.3 do not indicate the same 
difference. The RIA measurements were made as close as 
possible to the reflecting surfaces, 0.5 m above the ground 
plane and 0.2 m below the interaction plane, whereas in the 
PTA case the results represent an average of measurements 
made at puff positions at heights ranging between 1.0 and 
7.0 m. The results in Figures 8.2 and 8.3 therefore indicate 
that the difference in shock strength over the ground compared 
with that at the interaction plane may be dependent on the 
height above the ground - not an unexpected result. Determina- 
tion of Mach shock strength from measurements of the times of 
shock arrival at smoke puffs at various heights is difficult 
because an assumption must be made about the exact shape of the 
Mach shock front, in order to correctly assign shock radius 
values at smoke puff positions. At or near a reflecting 
surface the problem of shock shape is not so important as it 
is assumed that the shock is perpendicular to the surface. 
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Details of the problem in the PTA case and the manner in 
which the problem was dealt with for Shots 9 and 10 are 


described in Volume 1 of this report. 
3.4 Resolution of time histories 


The time histories of density and pressure shown 
in Figures 24 and 25 do not normally show a sharp rise at 
the shock front. This slow rise is not a real effect but 
| one inherent to the method of particle trajectory analysis, 
which does not permit a high resolution of density in space 
or in time because the average density of the air within a 
rectangular cell defined by four smoke puffs cannot be calcu- 
lated accurately until the shock has completely traversed the 
cell. The time of complete traversal may be as much as 5 ms. 
For the same reason the calculated time histories often 
anticipate the time of shock front arrival and do not resolve 
any weaker shocks subsequent to the first, although these 
shocks may be detected occasionally as a rounded bump in the 


normally exponentially decaying curve. Efforts are being 


made to improve the space and time resolution of density and 

7 pressure calculated from the particle trajectories. 

| The lack of resolution close to the shock front does not 
occur in the case of particle velocity, which can be measured 


with reasonable accuracy as soon as the shock has traversed 


the relatively small space represented by an individual smoke 
puff. This improved resolution is manifested also in the 
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dynamic pressure histories which depend on particle velocity 


squared. 


rs shbhi Srn 


3.5 Comparisons with Gauge Results 


The hydrostatic and total pressure time-histories were 
calculated at positions where gauges were located and the 
results from the particle trajectory analysis are compared 
with the corresponding transducer outputs in Figures 27.1 
€0Q\ 2754. 
| On the whole the agreement between the results from the 
two measurement methods is reasonable although as previously 
discussed the poor time resolution of the particle trajectory 
results does not permit indentification of multiple shocks. 
This is illustrated in Figure 27.1 for the 60:20 Location. 
However, the agreement between the two impulse curves is 
excellent. 

The agreement between the pressure curves at locations 
60:47, 60:50 and 60:53 in Figure 27.2 is not good. The cause 
of this has been traced back to the anomalous movement of a 
single smoke puff, at co-ordinates (1.9, 2.0) in Figure 13.2. 
It will be seen that this produces an unusually high density 
and pressure on one side of the puff and unusually low on the 
other side. Inspection of the film record confirms the 

anomalous movement of the single puff for which there is no 
Z ready explanation other than the possibility of a weak non- 


luminous jet from the lower charge (Patterson, et al., 1972). 
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The reasonable agreement between the total-pressure 
time-histories from the two measurement methods shown in 
| Figures 27.3 and 27.4 further confirms the validity of the 
technique used to calculate the total pressure, as described 
in the addendum to Volumes 1 and 2 of this report. 
In considering the above comparisons it must be remembered 
that determination of pressure time-histories was not an 
i objective of the particle trajectory analysis project, but the 
| reasonable agreement with gauge measurements gives some 
indication of the reliability of the method. Also the gauge 
measurements were made on the opposite side of the charges 
| to the smoke puff grid so that some differences might be 
expected due to slight non-symmetries of the blast waves. 
Although the effect is more clearly seen from the gauge 
Ei results rather than the particle trajectory analysis, it is 
interesting to note the difference in the blast wave signature 
at location 60:10, 10ft above the smooth ground, compared with 
that at 60:40, 10ft below the interaction plane. Above the 
ground the two initial shocks indicate that the gauge was 
above the triple point and observed the primary and the reflected 
shock. Below the interaction plane the single shock indicates 
that the gauge was in the Mach stem below the triple point. 
This confirms the photogrammetrical measurements (Dewey, et al., 
1975) which showed the triple point above the ground at this 
location to be at height of approximately 8.5 ft and to be 


approximately 10 ft below the interaction plane, i.e. coinciding 


with the gauge position. 


A more detailed comparison of the blast wave properties 
in the Mach reflection regions above the ground and below the 
interaction plane will be made in a subsequent volume of this 


report. 
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Introduction 


Volumes 1 and 2 of this report presented the results 
of the particle trajectory photogrammetry of Dipole West 
Shots 10 and 9, respectively. Those results included the 
variation with distance of the strengths of the primary 
spherical shocks from the pairs of charges and of the Mach 
shocks over the ground and ét the interaction plane between 
the pairs of charges; the particle velocity, density, 
hydrostatic overpressure and dynamic pressure fields throughout 
the blast waves at various times, and the time histories 
of these physical properties at several fixed positions. 
In particular, the time histories were given for positions 
which coincided with electronic gauge locations. At the 
time those reports were prepared the results of the gauge 
measurements were not available to the authors, and a 
direct comparison between the histories obtained from the 
particle trajectory analysis and those measured with electronic 
transducers was not possible. The gauge measurements have 
now been made available (Keefer & Reisler, 1975) and a 
comparison between the results obtained by the two measurement 
techniques for Shots 9 and 10 is the subject of this addendum. 
Similar comparisons for Shots 8 and 11 will be included in 


subsequent volumes of this report. 
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Two types of electronic pressure transducer were used 
in these experiments, arranged with the sensitive elements 
either side-on or face~-on to the blast, to measure the 
hydrostatic and total overpressures respectively. The 
calibrated signals from these transducers are presented by 
Keefer & Reisler (1975), together with the time history of the 
integral of each signal, representing the hydrostatic and 
total overpressure impulses. 

To facilitate the comparisons described above, the pressure-~ 
time histories at the gauge locations calculated from the 
particle trajectory analyses are plotted on identical scales 
to those used to report the gauge results. Computed hydro- 
static overpressures are plotted; also total pressures obtained 
by adding to the hydrostatic overpressures the computed 
dynamic pressures, with the application of a compressibility 
correction. Both types of result were also integrated to 
give the side-on and face-on pressure impulses. 

Comparisons were made for all gauge locations which lay 
within the smoke puff grids. The gauges were mounted on a 
vertical gun barrel at a distance of 60 ft (18.3 m) from GZ 


and at heights above the ground of 10, 15, 20, 27, 30, 33 and 


40 ft. (The interaction plane between the two charges in 


} Shots 9 and 10 was at a height of 30 ft.) 


Calculation of Total Pressure 


Dynamic pressures (%pu*, where p is the density and u 
y 


is the particle velocity) and hydrostatic overpressures 


| obtained from the analysis of the particle trajectories were 

used to compute total pressures after the application of a 

compressibility correction. The amount of the correction 

was computed as a function of local Mach number, the form of 

the function depending on whether the Mach number was 

greater than or less than unity. (If the local Mach number 
is greater than one, a bow shock forms around a pressure 
gauge and this further modifies the flow.) 

Two assumptions were made: the first was that the ratio 
of specific heats of air remains constant at 1.4. In other 
words, it was assumed that there were no real gas effects of 
importance. At positions on the 60ft gun barrel where the 
comparisons with the gauge results are made, the maximum total 
pressure computed was always less than 6 atmospheres. 


The second assumption was that the reflected pressure 


would not be detected because of the small size of the gauges 
used. The reflected pressure at a surface face-on to a blast 
wave lasts only until the pressure is relieved by the 
rarefraction wave produced as the shock defracts around the 
edge of the face-on surface. In the case of a small pressure 
transducer this relief of reflected pressure will be complete 


: in a few tenths of a millisecond. 
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For each point at which the total pressure was to be 
calculated the square of the local Mach number, M, was 


computed using 


where q is the dynamic pressure, S the absolute hydrostatic 
pressure and y the ratio of specific heats. The values of 

q and § obtained from the particle trajectory analysis are 
both ratios of the ambient atmospheric pressure. In the cases 
where M* was less than 1, the total pressure, T, was 


calculated using 


Y 
= 2 yaa 
a= 5 | Ooue +a 


In the cases where M? was greater than 1, the total-head 


pressure was calculated from 


It should be noted that the values obtained from the 
above equations are measures of absolute pressure. The 
plotted results are overpressure, namely T-l. The process 


described above was repeated for a sequence of times at each 


gauge location to provide the time histories. 
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Calculation of Pressure Impulses 


The hydrostatic and total pressure-time histories 
were integrated using the trapizoidal rule. Because of the 
lack of spacial resolution which is an inherent limitation of 
the particle trajectory analysis method imposed by the 
finite spacing of the smoke puffs, some of the pressure- 
time signals show a rounded leading edge. As discussed in 
Volume 2, this is an effect of the analysis method and not an 
indication of a real distortion of the pressure pulse. 
Therefore, the impulse integral at the leading edge was 
calculated by joining the peak pressure value (calculated 
from the shock velocity) to the maximum value in the pressure- 


time history obtained from the particle trajectory analysis. 
Results 


The results are presented in Figures 1 to 8. In each case 
the pressure and impulse curves obtained electronically 
(Keefer and Reisler, 1975) are shown as solid curves. The 
corresponding results obtained from the particle trajectory 
analysis are shown as dotted curves. For each such set of 
curves an arrow on the pressure axis indicates the peak value 


of pressure calculated from the shock velocity; that is, the 


initial value used in the impulse calculation. 


144 


Discussion 


The pressure and impulse curves obtained by the two 
completely different techniques, in general are in excellent 
agreement. The good agreement between the total pressure 
curves is particularly gratifying since this is the first time 
that such a comparison has been possible. This agreement would 
appear to validate both of the measurement techniques, and 
the compressibility corrections which were applied to the 
dynamic pressures. The hydrostatic overpressure results 

5 show agreement which is not as good and in this case one 
should suspect the particle trajectory analysis results as 
having the greatest error since hydrostatic pressure is the 
least accurate physical property obtained using the particle 
trajectory analysis technique. As explained in Volume 2, the 
finite spacing of the smoke puffs results in both a poor 

» definition of the pressure and density histories close to a 

shock front and the poor detection of weak subsequent shocks, 
which are clearly detected by the pressure transducers. These 
errors are not as great for the total pressure which is derived 
in part from the dynamic pressure. Although dynamic pressure 
involves the density, its major factor is the square of the 

‘ particle velocity, the most accurate measurement that can be 


obtained from the particle trajectory analysis. 
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Efforts are being made to improve the resolution of 
density and pressure obtained from the particle trajectory 
analysis. 

No attempt is made at this time to draw conclusions from 
the results presented here concerning the shock waves produced 
at the interaction plane and over the ground surfaces. Such 


matters will be considered in detail in a subsequent report. 
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The authors gratefully acknowledge the assistance of 
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